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Objective: The aim of this study was to investigate whether single nucleotide polymorphisms (SNPs) of dopamine receptor D2 
(DRD2) are associated with schizophrenia in Korean population. 

Metliods: Four SNPs (rs4648317, rs7131056, rs4936270, and rs1076562) of DRD2 were selected and genotyped by direct se- 
quencing in 197 schizophrenia patients and 370 control subjects. SNPAnalyzer, SNPStats, and Haploview version 4.2 programs 
were performed to analyze the genetic data. Multiple logistic regression models (codominantl, codominant2, dominant, recessive, 
overdominant, and log-additive) were used to evaluate the odds ratios (ORs), 95% confidence intervals (CIs), and p values. 
For multiple testing, p values {p^ were re-evaluated by Bonferroni's correction. 

Results: The genotype frequency of DRD2 rs4936270 SNP was associated with the development of schizophrenia (;C^0.0007, 
0R=1.71, 95% CI=1.16-2.52 in the codominantl model; p^Q.OW, OR=1.63, 95% 0=1.12-2.37 in the dominant model; /;=b.035, 
0R=1.41, 95% Cl=1.03-1.95 in the log-additive model). The allele frequency of rs4936270 was also associated with the develop- 
ment of schizophrenia (p=0.024, OR=1.45, 95% Cl=1.05-1.98). After Bonferroni's correction, the genotype distribution of rs4936270 
was still related to the development of schizophrenia (^"=0.0028 in the codominantl model; p''=0.044 in the dominant model). 
A linkage disequilibrium block consisted of rs4648317, rs7131056, and rs4936270. The CAT haplotype frequency was different 
between schizophrenia and controls [p^OOZQ). 

Conclusion: These results suggest that DRD2 SNPs may be associated with the development of schizophrenia in Korean 
population. 
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INTRODUCTION 

Schizophrenia is a complex psychiatric disorder char- 
acterized by a marked thinking disturbance and is also a 
devastating disorder with a high degree of morbidity.'' 
Twin, adoption, and family studies have implicated that 
schizophrenia is one of genetic disorders with high 
heritability.^'^' Several lines of evidences suggest that ge- 
netic factors contribute to the susceptibility of schizo- 
phrenia.'*'^' Recently, many genetic studies have been re- 
ported the relationship between schizophrenia and single 
nucleotide polymorphisms (SNPs) of candidate genes 
and/or copy number variants and susceptible loci through 
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the genome-wide association study.*"'^' 

The dopamine receptor D2 (DRD2) is a G pro- 
tein-coupled receptor inhibiting adenylyl cyclase activity. 
This receptor is located on postsynaptic dopaminergic 
neurons involved in reward-mediating mesocorticolimbic 
pathways.'' DRD2 consists of 2 isoforms, D2S (shot iso- 
form) and D2L (long isoform). D2L has an additional 29 
amino acids in the third intracellular loop, compared to 
D2S. D2L acts at postsynaptic sites and D2S plays a role 
in presynaptic autoreceptor functions. The pharmaco- 
logical effects of haloperidol are absent in D2L deficient 
mice, suggesting D2L is targeted by haloperidol.'"' The 
DRD2 signaling plays a crucial role in the physiopatho- 
logic conditions such as locomotion and drug abuse. It is 
well-known that DRD2 may be related to the suscepti- 
bility of schizophrenia. Antipsychotic drugs block post- 
synaptic DRD2. Postmortem studies have reported that 
DRD2 density is increased in the caudate, putamen, and 
nucleus accumbens of schizophrenia patients."' Drugs in- 
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creasing dopaminergic activity such as levodopa and bro- 
mocriptine aggravate schizophrenic symptoms. 
Although DRD2 is implicated in the pathogenesis of 
schizophrenia, the genetic determinants still remain. 

In this study, we investigated whether SNPs of the 
DRD2 gene are associated with schizophrenia in a Korean 
population. 

METHODS 



individuals assessed as mentally healthy through a general 
health examination program. Subjects with neurological 
diseases and any severe diseases were excluded. This 
study was carried out in accordance with the guidelines of 
the Helsinki Declaration and approved by the Ethics 
Review Committee of Medical Research Institute, Kyung 
Hee University Medical Center, Seoul, Republic of 
Korea. Written informed consent was obtained from each 
subject. 



Study Subjects 

Schizophrenia patients were enrolled among partic- 
ipants who visited at the departments of neuropsychiatry 
in the East- West Neomedical Center and Kyung Hee 
Medical Center (Seoul, Korea). A total of 197 Korean 
schizophrenia patients and 370 healthy subjects were 
recruited. The schizophrenia group (n=197) comprised 
113 males and 84 female, and the mean age of schizo- 
phrenia patients was 42.2±10.8 (mean±standard devia- 
tion) years. The control group (n=370) consisted of 187 
males and 183 females. The mean age of controls was 
44.3±6.3 years. Patients were diagnosed as schizophrenia 
by well-trained psychiatrists according to the Diagnostic 
and Statistical Manual of Mental Disorders (fourth edi- 
tion) using available information from interviews and 
clinical records. Control subjects were recruited among 
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Genotyping 

In this study, we focused on intronic SNPs (iSNPs) in 
the intron 1 region of DRD2: (1) we wanted to find out 
new SNPs related to schizophrenia susceptibility; (2) the 
methylation of CpG islands has been identified in the in- 
tron 1 of DRD2."^' We searched DRD2 SNPs in the NCBI 
websites (www.ensembl.org, www.broad.mit.edu/mpg/ 
tagger, www.hapmap.org, and www.ncbi.nlm.nih.gov/ 
SNP, BUILD135). SNPs with <5% minor allele fre- 
quency, < 10% heterozygosity, and monomorphic geno- 
type in Asian populations were excluded. Finally, we se- 
lected four iSNPs (rs4648317, rs7131056, rs4936270, and 
rsl076562) (Fig. 1). Peripheral blood samples were col- 
lected in EDTAor heparin tubes from each subject. DNA 
was isolated using the DNA Isolation Kit for Cells and 
Tissues (Roche, Indianapolis, IN, USA) and the geno- 
types of four SNPs were determined by direct sequencing 
(MACROGEN, Seoul, Korea). Polymerase chain re- 
actions (PCRs) using the sense and antisense primers of 
each SNP were conducted (Table 1). The PCR products 
were sequenced by an ABI PRISM 3730XL analyzer (PE 
Applied Biosystems, Foster City, CA, USA) and SeqManll 
software (DNASTAR, Madison, WI, USA) was used to 
analyze the sequencing data. 



Fig. 1. Gene map and location of single nucleotide polymorphisms 
(SNPs) In trie DRD2 gene. Exons ore marl<ed with box. The coding 
regions ore blacl< boxes and untronsiotion regions ore white 
boxes. Arrow indicates the location of each SNP. 



Statistical Analysis 

SNPStats (ht^://bioinfo.iconcologia.net/index.php?module 
=Snpstats) and SNPAnalyzer Pro (ISTECH, Goyang, 



Table 1. Primer sequences of DRD2 SNPs 



, SNP 



Sequence (5'-3') 



DRD2 



rs4648317 

intron 
rs 1076562 

intron 
rs4936270 

intron 
rs71 31056 

Intron 



Forward 
Reverse 
Forward 
Reverse 
Forward 
Reverse 
Forward 
Reverse 



GGTACAGGGGATTTCCAAGGAG 

AGAGGCACATGGACACTGCTGG 

TAHGCAACCTCCCCAGACTCA 

CAACAGCCAAHAATCAAGCAA 

TCHAGGACTGGTGGGAGGAAG 

AAGAAAGTGCTGCFGHCCAG 

GCTTTCTTGCCCACAGTACATT 

GCAGGGGGAAATAAGAGGCAAG 



306 
391 
360 
361 



DRD2, dopamine receptor D2; SNP, singe nucleotide polymorphism; bp, base pair. 
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Korea) were used to calculate odds ratios (ORs), 95% con- 
fidence intervals (CIs), and p values adjusting for age and 
gender as covariables. Multiple logistic regression models 
were employed to analyze genetic data: codominantl 
(major allele homozygotes vs. heterozygotes), codomi- 
nant2 (major allele homozygotes vs. minor allele homo- 
zygotes), dominant (major allele homozygotes vs. hetero- 
zygotes+minor allele homozygotes), recessive (major al- 
lele homozygotes+heterozygotes vs. minor allele homo- 
zygotes), overdominant (major allele homozygotes+mi- 
nor allele homozygotes vs. heterozygotes), and log-addi- 
tive (major allele homozygotes vs. heterozygotes vs. mi- 
nor allele homozygotes). The linkage disequilibrium (LD) 
block and haplotypes were evaluated using Haploview 4.2 
(Daly Lab, Cambridge, MA, USA). Bonferroni's correc- 
tion was applied by multiplying the p values by the num- 
ber of SNPs. Hardy- Weinberg equilibrium (HWE) was 
calculated using the chi-squared test. Statistical sig- 
nificance was evaluated using SPSS ver. 15.0 (SPSS Inc., 
Chicago, IL, USA) and set at;? < 0.05. 



RESULTS 

Genotype Analysis 

The genotype and allele Irequencies of the four exam- 
ined SNPs are shown in Table 2. Four iSNPs (rs4648317, 
rs7131056, rs4936270, and rsl076562) of DRD2 were in 
HWE in the schizophrenia and control groups, respec- 
tively (p > 0.05, data not shown). Three SNPs (rs46483 1 7, 
rs7131056, and rsl076562) were not associated with 
schizophrenia. However, the genotype frequency of 
rs4936270 was associated with the development of schiz- 
ophrenia (p=0.0007, OR=1.71, 95% 01=1.16-2.52 in the 
codominantl model;/?=0.011,OR=1.63, 95% 01=1.12-2.37 
in the dominant model; jc=0.035, OR=1.41, 95% 01= 
1 .03- 1 .95 in the log-additive model). The allele frequency 
of rs4936270 was also associated with the development of 
schizophrenia (p=0.024, OR=1.45, 95% CI=1.05-1.98). 
The C/T genotype and C allele Irequencies of rs4936270 
was higher in the schizophrenia group (34.7 and 20.3%) 
than in the control group (23.5 and 15.0%). After 
Bonferroni's correction, the genotype distribution of 



Table 2. Genotype and allele frequencies of DRD2 SNPs In control subjects and sctnizoptnrenia patients 



SNP 




Type 


Control 




SPR 


Model 


OR (95% CI) 


P 


Bonferroni' 












correction 


rs4648317 


Genotype 


C/C 


133 


36.0 


71 


36.2 


Codominantl 


0.99 (0.67-1.46) 


0.96 


1.00 


Intron 




C/T 


175 


47.3 


95 


48.5 


Codominant2 


0.91 (0.54-1.55) 


0.74 


1.00 






T/T 


62 


16.8 


30 


15.3 


Dominant 

Recessive 
Log-additive 


0.97 (0.67-1.40) 
0.92 (0.57-1.48) 
0.96 (0.75-1.24) 


0.87 
0.73 
0.77 


1.00 
1.00 
1.00 




Allele 


C 


441 


59.6 


237 


60.2 




1 










T 


299 


40.4 


157 


39.8 




0.98 (0.76-1.25) 


0.86 


1.00 


rs 1076562 


Genotype 


A/A 


95 


25.7 


54 


27.6 


Codominantl 


0.96 (0.63-1.46) 


0.85 


1.00 


Intron 




A/G 


181 


48.9 


100 


51.0 


Codominant2 


0.82 (0.50-1.35) 


0.44 


1.00 






G/G 


94 


25.4 


42 


21.4 


Dominant 
Recessive 
Log-additive 


0.91 (0.62-1.36) 
0.84 (0.55-1.28) 
0.91 (0.71-1.16) 


0.65 
0.42 
0.45 


1.00 
1.00 
1.00 




Allele 


A 


371 


50.1 


210 


53.3 




1 










G 


369 


49.9 


184 


46.7 




0.88 (0.69-1.13) 


0.31 


1.00 


rs4936270 


Genotype 


C/C 


271 


73.2 


122 


62.2 


Codominantl 


1.71 (1.16-2.52) 


0.0007 


0.0028 


Intron 




C/T 


87 


23.5 


68 


34.7 


Codominant2 


1.05 (0.38-2.89) 


0.93 


1 






T/T 


12 


3.2 


6 


3.1 


Dominant 
Recessive 
Log-additive 


1.63 (1.12-2.37) 
0.89 (0.32-2.43) 
1.41 (1.03-1.95) 


0.011 
0.82 
0.035 


0.044 

1.00 
0.14 




Allele 


C 


629 


85.0 


314 


79.7 




1 










T 


111 


15.0 


80 


20.3 




1.45 (1.05-1.98) 


0.024 


0.96 


rs71 31056 


Genotype 


C/C 


121 


32.7 


62 


31.6 


Codominantl 


0.98 (0.65-1.46) 


0.91 


1.00 


Intron 




C/A 


174 


47.0 


92 


46.9 


Codominant2 


1.09 (0.66-1.77) 


0.75 


1.00 






A/A 


75 


20.3 


42 


21.4 


Dominant 
Recessive 
Log-additive 


1.01 (0.69-1.47) 
1.10 (0.72-1.69) 
1.03 (0.81-1.32) 


0.96 
0.66 
0.78 


1.00 
1.00 
1.00 




Allele 


C 


416 


56.2 


218 


55.3 




1 










A 


324 


43.8 


176 


44.7 




1.04 (0.81-1.33) 


0.78 


1.00 



DRD2, dopamine receptor D2; SNP, singe nucleotide polymorphism; SPR, schizophrenia; OR, odds ratio; CI, confidence interval. 
The p values were calculated from logistic regression analysis adjusting age and sex. Bold number indicates significant association. 
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rs4936270 showed statistical difference between schizo- 
phrenia and controls (p'^=0.0028 in the codominantl mod- 
el; p'^=0.044 in the dominant model). These results suggest 
that rs4936270 may be a risk factor for schizophrenia 
development. 

Linkage Disequilibrium and Haplotypes 

In order to evaluate the LD block and haplotypes of the 
four tested SNPs, we used Haploview 4.2 program. In Fig. 
2, a LD block was made among rs4648317, rs7131056, 
and rs4936270 SNPs in the control group (rs46483 17 and 
rs7131056, D'=l, r^=0.528; rs7131056 and rs4936270, 
D'=l, r^=0.227; rs4648317 and rs4936270, D'=l, 
r^=0.120). We also investigated the association between 
the haplotypes and schizophrenia. As shown in Table 3, 
frequency in four haplotypes showed >0.1. The dis- 
tribution of TCC, CAC, CAT, CCC haplotypes was 0.400, 
0.275, 0.166, and 0.157, respectively. The CAThaplotype 
was weakly associated with the development of schizo- 
phrenia (p=0.039). In addition, the genotype frequencies 
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Fig. 2. Linkage disequilibrium block consists of rs4648317, rs7131056, 
and rs4936270. 



of the four examined SNPs (rs4648317, rs 1076562, 
rs4936270, and rs7131056) were compared to those of 
European and other East Asian populations. The genotype 
frequencies of three SNPs (rs4648317, rs 1076562, and 
rs713 1056) were similar to those of Japanese, whereas the 
genotype frequencies of rs4936270 were similar to those 
of European (Table 4). 

Sample Power 

The power of the sample size was estimated using a ge- 
netic power calculator (http://pngu.mgh.harvard.edu/~purceU/ 
gpc). When the conditions of the sample power (a =0.05, 
risk=two-fold, the number of case=80% power) were em- 
ployed, we had 0.9 143 forrs4648317 (n=138), 0.9251 for 
rs7131056 (n=133), 0.9184 for rs4936270 (n=136), and 
0.8906 for rsl076562 (n=151). Therefore, our case-con- 
trol study was sufficiently powerful to study a positive 
correlation. 

DISCUSSION 

A number of studies have been pubhshed the relation- 
ship between DRD2 SNPs and schizophrenia. Betcheva et 
al."^ demonstrated that rs6277 (C957T, Pro319Pro) was 
significantly associated with schizophrenia after the repli- 
cation study using independent sample set (OR=1.76, 
p=0.001). Fan et al}^^ reported that T allele of rs6277 
(Pro319Pro) was weakly associated with schizophrenia 
(0R=1 .58, 95% CI=1 .03-2.43,^7=0.034), while the geno- 
type and allele frequencies of rs6275 (C939T, His3 1 3His) 
were not associated with schizophrenia in Chinese Han 
population. They suggest that rs6277 T allele may play a 
role in the genetic vulnerability for schizophrenia. Amis- 
sense SNP rsl801028 (Ser311Cys) was reported to have 
an effect on risk for schizophrenia. In conttast, Glatt and 
co-workers^"' did not find a significant effect of rs 1 80 1 028 
for association of schizophrenia in a family-based study. A 
study using in Southern Indian population found an asso- 
ciation between the genotype disfribution of rs6275 and 



Table 3. Haplotype analysis of DRD2 SNPs in control subjects and sctiizoptirenia patients 



Control SPR 

Haplotype Frequency Ctii square p 

TCC 0.400 298.6 441.4 154.7 239.3 0.13 0.72 

CAC 0.275 213.4 526.6 98.3 295.7 1.95 0.16 

CAT 0.166 110.6 629.4 77.7 316.3 4.23 0.039 

CCC 0.157 117.0 623.0 61.0 333.0 0.02 0.88 



DRD2, dopamine receptor D2; SNP, singe nucleotide polymorptiism; SPR, schlzoptirenia. 

Haplotypes comprised of rs4648317, rs7131056, and rs4936270. Bold number Indicates significant association. 
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Table 4. Genotype and allele frequencies of DRD2 SNPs in European and Asian populations 



SNP 




Type 


Korean* 


CEU (%) 


iHCB (%) 


JPT (%) 






n 


% ftM 




rs4648317 


Genotype 


C/C 


133 


36.0 


75.2 


41.9 


44.2 




intron 




C/T 


175 


47.3 


23.0 


37.2 


40.7 








T/T 


62 


16.8 


0.02 


20.9 


15.1 






Aileie 


C 


441 


59.6 


86.7 


60.5 


64.5 








T 


299 


40.4 


13.3 


39.5 


35.5 




rs 1076562 


Genotype 


A/A 


95 


25.7 


5.3 


16.3 


26.7 




Intron 




A/G 


181 


48.9 


36.3 


58.1 


50.0 








G/G 


94 


25.4 


58.4 


25.6 


23.3 






Allele 


A 


371 


51.2 


23.5 


45.3 


51,7 








G 


369 


49.9 


76.5 


54.7 


48.3 




rs4936270 


Genotype 


C/C 


271 


73.2 


76.1 


55.8 


53.5 




intron 




C/T 
T/T 


87 
12 


23.5 
3.2 


23.9 


41.9 
2.3 


38.4 
8.1 






Aiieie 


C 


629 


85.0 


88.1 


76.7 


72.7 








T 


111 


15.0 


11.9 


23.3 


27.3 




rs7131056 


Genotype 


C/C 


121 


32.7 


18.6 


18.6 


27.9 




Intron 




C/A 


174 


47.0 


48.7 


53.5 


48.8 








A/A 


75 


20.3 


32.7 


27.9 


23.3 






Aileie 


C 


416 


56.2 


42.9 


45.3 


52,3 








A 


324 


43.8 


57.1 


54.7 


47.7 





DRD2, dopamine receptor D2; SNP, singe nucleotide polymorphism; CEU, Utah residents with Northern and Western European ancestry 
from the CEPH collection; CCB, Han Chinese In Beijing, China; JPT, Japanese in Tol<yo, Japan. 
*Controi subjects in this study (http://www.ncbi.nlm.nlh.gov/projects/SNP). 



schizophrenia (x^=8 .9 1 , jc=0.0 1 1 ).^'' The rs2734839 iSNP 
was reported to have an association with schizophrenia 
with late onset age.^^^ Lee et alP'' showed that rsl799978 
(promoter SNP), rs4245149 (iSNP), and rs4620755 
(iSNP) were not associated with schizophrenia in Korean 
population. 

Focusing on the our tested SNPs, rs46483 17 was asso- 
ciated with better performance on the stop task in the drug 
free condition and lower scores on the impulsivity 
subscale.^'*' Laucht et al.^^^ found that carriers of the T al- 
lele of rs46483 1 7 were higher rates of current smoking 
and scored higher on nicotine dependence than their al- 
lelic counterparts. Volter et al.^^^ reported that rs46483 17 
was related to prepulse inhibition of the acoustic startle re- 
sponse in the Munich sample, but not in the London sample. 
A Finnish population-based cohort showed an association 
between rs7 13 1056 and social phobia 0=0.0084).^'' Todt 
et al^^^ reported that rs7131056 was associated with mi- 
graine with aura in a large German migraine with aura 
case-control sample (uncorrected allelic />=0.0018, 
OR=1.28). However, this correlation did not conform in 
2,937 British control individuals. The rs7131056 iSNP 
was reported to have initial significant association with 
nicotine dependence in European Americans, but not after 
correction for multiple testing, indicating a weak associa- 



tion of rs7131056 with nicotine dependence. Xing et 
al.^"'' reported that rsl076562 did not show any associa- 
tion with the response to risperidone in Chinese patients 
with schizophrenia. 

A famous SNP as known as TaqI DRD2 (rs 1 800497) is 
located within exon 8 of the adjacent gene, ankyrin repeat 
and kinase domain containing 1 (ANKKl). Therefore, we 
did not discuss about rs 1800497 studies in this study. 

In conclusion, we investigated whether four iSNPs 
(rs4648317, rs7131056, rs4936270, and rsl076562) of 
DRD2 are associated with schizophrenia in Korean 
population. We found that the genotype distribution of 
rs4936270 was associated with schizophrenia (p''=0.0028 
in the codominantl model; in the dominant mod- 

el) and the CAT haplotype was different between schizo- 
phrenia and controls (p=0.039). These results suggest that 
rs4936270 of DRD2 may be associated with the suscepti- 
bility of schizophrenia. To confirm our results, replication 
studies of a large number of cases and/or additional stud- 
ies of other DRD2 SNPs will be required. 
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